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ABSTRACT 


Velocity and temperature fluctuation measurements were 
Made over the open ocean from instruments mounted on the 
R/V ACANIA. These data were examined to determine the 
validity of Br eset formulations and prediction techniques. 

Values of momentum flux, cia: were inferred from the 
rate of dissipation of turbulent kinetic energy, €. Dissi- 
pation values were obtained from spectral estimates and 
Pinte rsca lS ects ti Mates m savallcs of su were examined for 
representativenessS on the basis of the constant-flux assump- 
tion and by comparisons with other studies. The vertical 
variation of the dissipation rate was examined for possible 
effects of stability and wind-wave coupling. 

The momentum flux, computed from spectra, supported the 
constant-flux assumption for neutral conditions. For 
Der moOdSmoremanstablility, —themstability corrections jappli ed 
Eo enemverticlal viarlation of € resulted in the’ proper VWad=- 
Wastmentse toward the predicted slope. “I[he reductions in 
momMmemtcum transfer during periods GEystable stratifitecations 
were consistent with wind-wave coupling effects described 
ieee avidson. The shape of a Spectrum Of temperature Eluc-— 


tuations was in agreement with predictions. 
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Le. INTRODUCTION 


The structure of turbulence in the marine boundary 
layer has become of increasing importance in naval systems 
design. To gain a better understanding of turbulence sta- 
tistics, such as the spectral distribution of wind fluc- 
tuations, the effects of wind-wave coupling and the 
eau choy of present prediction expressions, personnel at 
the Naval Postgraduate School, Monterey, California have 
undertaken a study to obtain the shipboard observational 
data necessary Foe accurate examination of these prob- 
lems. The group consists of personnel from several depart- 
ments; Meteorology, Mechanical Engineering, Oceanography 
and Physics, whose primary objective is to obtain observa- 
tional data and then to evaluate and improve existing pre- 
diction techniques. 

Several observational experiments have been conducted 
andmwere followed by data Oe xass and analyses periods. 
The second observational experiment is the one considered 
Mr chas Study and wai Conducted from 18 September 1973 to 
ec ocpuember 1973. Meteorological measurements were made 
from the oceanographic research ship, the R/V ACANIA shown 
Miererqurce loan the vicinity of San Nicolas Island. Fig- 
unre 2 illustrates the “near" oc@an environment and the 
ACANIA'S Anchorage where the majority of measurements were 


ak Cone, 








Figure 1. The R/V ACANIA is operated by the Department 
of Oceanography, NPS. 
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R/V ACANIA 


San Nicolas 


Island 


Figure 2. San Nicolas Island and ACANIA's Anchorage. 
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The meteorological measurements were of three differ- 
ent types. Two of theSe pertained to different scales of 
measurements and the third was of the influence of the 
moving platform on mean and turbulent measurements. The 
scales considered were the turbulent scale and that which 
is associated with the mean wind and thermal stratifica- 
tc 1TOr. 

The scope of this study is confined to turbulent fluc- 
tuations of velocity and temperature, although reference 
will be made to mean wind and thermal stratification mea- 
surements. The latter are important since they are re- 
quired to relate to results from other experiments or for 
extrapolating to regions or conditions where turbulence 
measurements are not made. A nearly unique problem facing 
the Navy is that present expressions for the turbulent 
processes are empirical in nature and have only been pre- 
viously validated, or at least tested, overland. 

The purpose of this investigation is to test several 
expressions in the "near" ocean environment, interpret the 
results and evaluate existing turbulence expressions. To 
aceoumptlamshi sthisS, lt was adecided, tosconcentrate thegstudy 
On Specitiempormu iiationswee tnemrormuilataons ineiude 4i)) an 
evalmatmon Of the conSitteant Elux, Lalyver (Zz) “the variations 
Sire ene VISCOUS GdiSSipation ct turbulent Kinetic energy) 
with height and, where data permitted, (3) the variation 


of temperature fluctuations. 
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Ii. OBSERVATIONAL AND THEORETICAL CONSIDERATIONS 


A. GENERAL BACKGROUND 

This study represents one of the first attempts to 
Make simultaneous, multiple-level ives Pa: over water 
from a ship. Previous measurements have been primarily 
overland or over water from near stable platforms. A 
study which closely approximates this experiment was de- 
scribed by en a et al Gy a) The measurements re- 
ported in that study were made from the Scripps Floating 
Instrument Platform (FLIP) during the Barbados Oceano- 
graphic and Meteorological Experiment (BOMEX). From FLIP, 
velocity and temperature measurementsS were made with a 
roving probe which cycled continuously between measurement 
Stations located 2, 4, 7 and 12 meters above mean Sea 
level. Results in the present study will be compared to 
those obtained by Stegen, et al, as well as to other con- 
temporary studies. 

The following paragraphs will provide a review of the 
theoretical considerations needed to interpret the data. 
Paoms thiesie conisideraitions, the fraction vielochty (u,) will 
be estimated from values of the dissipation of turbulent 
minetremenergy (6) and € willl be, obtained Trom "spectra: 
For this reason we must examine both the turbulent kinetic 


energy balance expression and Kolmogorov's hypotheses, 


igs 





B. TURBULENT BALANCE EXPRESSIONS 
A general expression for the turbulent kinetic energy 


balance has the form 




















a vs du g Pp aS 1 
ye gaa ox tai 73° ime rox, 895 eo —)! 
VV 
d = Thay ea 


If all eddy transports except those along the vertical are 
neglected, the turbulent kinetic energy balance can be ex- 


pressed as 














ew, 


¢) = Ul ————— 
Dt | 2 debe one OX. 


(2) 
The time variation of turbulent kinetic energy is thus Seen 
to depend on four processes represented by the four terms 
On §thef right hand shde of Eiguiation (2). The finmstsrepre= 
sents production of turbulent kinetic energy by the inter- 
action of the Reynolds stress with the mean shear flow. The 
second represents turbulent kinetic energy released or 
gained due to the work of vertical motions against the 
effects of buoyancy. The third represents viscous dissipa- 
tion of turbulent kinetic energy. The final term represents 
bomieemessube wWOrk and the £luxPdivergence of turbulent 


kKimteure  wenergy due to tramsfer by turbulent motion. 
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To obtain those expressions normally used for estimat- 
. —— 2 
ing the shear stress or eddy flux (-vViV, =30, 7eor the 
basis of turbulent kinetic energy balance, the following 


assumptions and definitions are made: 








1) assume that the turbulent kinetic energy is in 
a statistically steady state 
: 2 
V 
) i = 
ar 5 ) = 0 
ii) conditions are horizontally homogeneous for the 


irregularity fields (1.e., spatial variations 
of mean fluctuations are taken to be negligible 
over distances comparable to the largest scales 
of interest in the analysis). 


0 O 
oe 
1 


) = ox, | ) = 0 


et 1) conditions are near neutral 


G ~ 
=“— v.90 = 
T 3 : 
O 
iv) define the coordinate system so 
u =u = O 
ey ana 


The above assumptions and definitions lead to the follow- 


ing turbulent energy balance expression 





= 
© 
—s 
fe 
NO 
OJ 
[od 
a 


(3) 





Furthermore, if the divergence terms are negligible, 


Equation (3) becomes 








SY = -€ (4) 
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which states that mechanical production of turbulent 
kinetic energy equals dissipation of turbulent kinetic 
energy. 


Equation (4) can be rewritten as 


2 9u 
u — = € 
3 
2 <a ’ * . 
where u, = “ViV¥, a For neutral conditions in a constant- 


stress boundary layer the following is an accepted expres- 


Samon for the wind’shear 


€ 


Rei gy (6) 


where k is Von Karman's constant and z is the height. 
Combining Equations (5) and (6) yields the following esti- 


iitomwoG Me DbDascad son Measumenents  OLec eat Known srelghie a2. 


u, = uae eo 


* 


(7) 


Equation (7) can be used to both estimate the value of the 
momentum flux, Ck ae and to examine if it is constant with 
respect to height, the constant~-flux assumption, 

Tablie I provides) the momentumyflux calculations; pre- 
sented by Stegen, et al with € values derived from time 
derivativels of the velocity fluctuations. These values 
will be compared to those obtained from shipboard measure- 
ments to establish the validity of the expressions and 


assumptions used in this study. 
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TABLE I 


Results of Momentum-Flux CellC omiaons 


Height kze OUR u, (corrected) 
(m) (cm3sec73) (cm sec™1) (cm sect) 


3 


el. Oo 4.02x10 oe 14.0 
6.70 ero or on 1S ural, 4 age 
Syegehe 6.36x10> Loe 17.4 
APE ee 6.03x10° 18.2 17.4 


TResults taken from Stegen, et al (1973) 


ey, 





Equation (7) can be rewritten as 
S 
De ee wee edie Jk (8) 


Therefore a plot of Ln z versus Ln € should yield a 
straight line with -1l slope. The intercept of this line 
with the Ln Zz axis, 2 = l, will yield the estimate of 


u A plot of this shown by Stegen, et al appears in 


*° 
Figure 3. Deviations from this particular relation could 
be interpreted for possible deviation from the constant- 


flux assumption or for influence of wind-wave coupling 


Or non-neutral stratifications. 


Cz SIMILARITY EXPRESSIONS 

The parameter € can be estimated on the basis of 
Similarity expressions for the inertial subrange. Accord- 
ing to Kolmogorov's second hypothesis, the inertial sub- 
range is a region (in wave number space) where there is 
negligible dissipation of energy and where inertial energy 
transfer to higher wave numbers is the dominate feature, 

hagutem@4) From Lumley vanad = Panotsky (1964) depicts the 
inertial subrange considering spectral energy transfers. 

Based on a dimensional analysis Kolmogorov's second 


nypotnesmse yields the following functional relation 


E(k) = eid = 


where E(k) represents the three-dimensional energy density, 


k aS the wave number and ¢c iS an empirical constant. 
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o 18M 1130 
© CDC 3600(1) 
ACDC 3600(2) 


@ Ave Value 





Figure 3. Measured Values of € vs z, the Height Above 
the) Mean #Sea@sSuriace (alften Stegen, etal). 
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Gain E(k) 


Inertial 
Transfer 





Viscous 
| Dissipation 


Loss ——— 


Figure 4. A Schematic Drawing Showing Spectral Transfer 
sandeillustrateng AF) Energy Containing Region 
(B) Inertial Subrange (C) DisSipation Range 
and (D) Universal Equilibrium Range. 
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Unfortunately, E(k) is impossible to measure. It can be 
Shown, however, that the one-dimensional energy density, 


91, (*), which can be measured, also has the form 


9 2/3,7973 


(kK) = a (9) 


ra 


where @ = constant. 

Furthermore, measurements at the scales of interest 
have to be made at a fixed point in the flow and the spec- 
tra therefore are obtained for temporal frequency, f. In 
order to use Equation (9), time and space scales must be 
related using Taylor's (1938) "frozen turbulence" hypo- 
Ene sis: 


k = 21f/U 


where U is the mean wind speed, measured by a cup anemo- 
meter. The term "frozen turbulence’ implies that the 
Bur bulence Sremains unchangededuring the time required for 
it to sweep past the probe. 

mires tol lowing “form. of -Equmition )(9) 1S often utwlized 


and will be here, 


27,3273 
52/3 =2/ 


ae Me = psd, 0) ae (10) 


11 


Them parameter © willl be obtained fiom the spectra using 
Equation (10). The € values obtained from simultancous, 
multiple level measurements will then be used to estimate 


2 
Eweomvaluceot sthe monentunm flux, ul, 9, from Equiqtion (7). 


pat 





There is another way of estimating € based on the 
scale associated with the upper limit of the inertial sub- 
range. If the spectrum is in statistical equilibrium it 
is independent of the source of the turbulence and depends 
only on €, Y and ke. A dimensional analysis yields the 


following expression for the inner scale, ly 


TOM=BiCy eye 


A (11) 


where kG 271 /u) is the dissipation wave number and y 

is the kinematic viscosity of air. The dissipation wave 
number can be obtained from spectra which show scale 
separation between pure inertial transfer (-5/3 regions) 
and the viscous dissipation region. Dissipation of turbu- 
lent kinetic energy can then be obtained from Equation (11) 


and compared to those obtained from Equation (10) to estab- 


lish greater confidence in the spectral estimate. 


D. MURDER THEORETICAL CONSIDERATIONS 

Although the expressions previously described are based 
on accepted and consistent theories, there are several 
serious assumptions which must be considered in applying 
them to observational data. In many instances the theory 
is empirical and the formulating analyseS were based on 
Sve hangdedatan  Aarecent sefftort in turbulence theory has 
be CummbtompYrovesOrsdisprovemthe validity of overland formu- 
lations for non-neutral conditions or over water regimes 
where the wave influence may be important. The latter has 


beenmaescribed by Davidson (1974). 


eZ 





— 


In developing the turbulent kinetic energy balance, 
given by Equation (3), the existence of meutral conditions 
waS a critical assumption. For non=-neutral stratifica- 
tions, Equation (7) would include a buoyamcy term and be 


of the form 


3 
ar Za Zi 
Ee = — _ = 
ka [O, (>) *) (12) 
ona 
3 
uy | 
= i et 
mas 3 aon . — 
Wie Wem hee TiS /gkV 8 = Monin~Obuhkov Bength; V8 = heat 
flux, 
Ri = Richardson Number, 
Py =Sempirical ##unction of i r 
?, = empirical function of Ri. 


Potaeete  non-neutral Wcasée jetie expression corresponding to 
Equation (8) becomes 


3 
uy, 





Ln € = -Lm z@+ Ln @, (Ri) + Ln (13) 


Possible deviations from the constant flux assumption or 
wind wave coupling influence exhibited by spectral results 
could be examined by testing overland expressions for 9. (Ri) » 
PSG pimleynCOrrectlonecOUlLdsthen wee applijred for periods of 
non-meutial stratifications. A possible effect of the 
Pistaweleyeonmsthe predicted plots of Ln E€ versus Ln z is 


Pelsthaweam@ain Figure. 5.. 
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Ln( Ux/ k) 


oO 







neutral case 
iné 3 
e£= Ux/kZ 
in Z 
3 
Ln( Ux/,) 
in€ non-neufral case 


correction 


ae Ln ). (Ri) 


Ln Zz 


Figue@emow ean Lllwstmacion of a Possible Effect of 
Saab ieiteyewom CllemvVariatzon of € with Neight. 
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FONG OBSERVATIONS 


A. THE EXPERIMENT 

The data considered in this study were obtained from 
measurements made aboard the R/V ACANIA, which is operated 
by the Department of Oceanography of the United States 
Naval Postgraduate School. AS previously stated, San 
Nicolas Isiand was the experimental site and the majority 
of data were collected while te Ship was at anchorage near 
the northwest tip of the island. Additional data were col- 
lected while the ship was enroute between Monterey and San 
Nicolas Island and also while the ship was transiting to 
the west side of the island. Since the datameo Pilec ted =:rom 
anchorage comprise the largest percentage of successful 
runsjand nepresents all stability conditions, data runs 
selected for analysis in this study are primarily anchorage 
datacme tablet li@iwse a listing of the various data runs’) ex- 
amined. Figure 6 shows snow simi 4 traick jin themvicimiity of 
San Nicolas Island and Figure 7 shows the route to and from 


the island. 


B. DE onl let ON oie EOULP MENG NDSMEASUREMENTS 
Figure 8 shows the sensor placement at various levels 
along with the variables measured at each level. 
Mee velocity Measurements 
Velocity fluctuation measurements were made at three 


levels on the forward mast and one level on a Single bow 
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Data Periods Considered in This Study 


Date Time Digital Channel Parameter 


20 Sept. 73 1103-1138 
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1 ohe subscript refers to the level at which the 
measurement was made. Level 1 is the uppermost 
level. 
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TABLE II (continued) 
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Date Time Digital Channel Parameter 


20 Sept. 73 2323-2341 
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Figure 7. The San Nicolas Island Region. 
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mast as shown in Figure 8. Thermo-Systems (TSI) Model 
1210 probes with sliding support shields were selected 
for these measurements. The sensor consisted of a plati- 
num coated tungsten wire 0.00015 inches in diameter and 
0.060 inches long, small enough to resolve the viscous 
dissipation scale without needing to apply wire length 
corrections. The wire was fitted with plated ends for 
isolating the sensing area and thereby minimizing flow 
disturbance. 

Electronics associated with Pe: probe were a TSI 
Model 1054B Linearized Anemometer, a TSI Model 1056 Variable 
Decade Module and a Model 1057 Signal Conditioner. Schema- 
tics of the anemometer circuits appear in Appendix A. The 
anemometer had a linear frequency response from DC to 10 KHz 
and the variable decade module operated with a 0-60 ohm 
range. The signal conditioner had a linear frequency 
response from DC to 400 KHz and permitted voltage suppres- 
STOnoLn one volt steps from,0—29 volts. Tie voltage outpuc 
from the conditioner was SORE on magnetic tape. 

Sensor placement required 200 foot cables areca 
Gtethe usual 15 foot cables. tThis wesulted jan a decrease 
ims y st ememrequency response butmehe sdecrease had little 
ettect, in the Erequency band Of intlerest. Ten probes were 
avallables for the Se and they were placed into 
BOUT GrEOUDS oa Each group Gf probes was calibrated with 
separate cables and electronics which permitted a complete 


Svisumenmcalibration fLorveach of the wfour ‘semsing levels. 


ey 





The calibration was completed with a manometer and a TSI 
Model 1125 Calibrator shown in Figures 9 and 10. Empiri- 
cal relations used in interpreting the hot wire signals 
are presented in Appendix B along with the computer code 
used to evaluate these measurements and resultant calibra- 
tion curves for each probe. 

The performance of the hot-wire system was excel- 
lent throughout the experiment. Only one probe failure 
was experienced despite continued probe exposure during 
ther d@ta’tgathering period. 

Ze Temperature Measurements 

Temperature fluctuation measurements were made at 
two levels on the forward mast of R/V ACANIA as Shown in 
Figure 8. Electronics used were TSI Model 1044 Tempera- 
ture Modules. The unit power supply consisted of six volt 
batteries. It can be utilized to make both temperature 
and AT measurements. The temperature module is designed 
for measurements of temperature variations down to 0.001°C 
qc. fLrequienciles up tomel) KHz. The schematics for the Model 
1044 Temperature Module are shown in Appendix C. 

Prebe £Lailurme occurnved.at the beginning of ithe 
first data run and to permit temperature measurements to 
Smeanwesthie probe unit wasslittedgwith a 60 ohm platinum 
wire sensor. 

3. Parallel Measurements 
Mean measurements of the wind speed, temperature 


and relative humidity have been described by Cavanaugh 
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(1974). Measurements of ship's motion have been described 


by Welsh (1974). 


Cx SIGNAL PROCESSING AND RECORDING 
1. Analog Recording 

All data were recorded simultaneously in analog 
form with a 14 channel Sangamo Model 3562 FM Recorder. 
All recordings were made at a recording speed of 7 1/2 
inches per second. At this recording speed the recorder 
had a frequency range from DC to 25 KHz. 

A schematic diagram of the complete recording sys- 
tem is shown in Appendix D. The amplifier and suppressor 
combination in series with the sensor electronics and re- 
corder channels were necessary to improve the signal to 
noise ratio of the inputs. The hot-wire signals were high 
pass filtered at 2 Hz during the majority of the runs. 
Two data runs were made without the filter to permit low 
frequency correlation between the velocity signal and 
ship's motion. The results of this correlation are dis-~ 
cussed by Welsh (1974). The temperature signals were low 
pass @Giltieredeat 500 Hz. In addition to thie signal re- 
cordings, a voice track was utilized to record observa- 
Edonse pertinent to the @ata runs. Constant monitoring of 
amplified and preamplified signals was performed with an 
oscilloscope to evaluate the signals and to insure that 
the signal would not exceed amplifier or recorder input 


limits. 


a0 


2. Real Time Analysis 

Real time analysis was accomplished uSing a Spec- 
tral Dynamics Corporation (SD) Model SD301C Real Time 
Analyzer. This ee was used in conjunction with a Model 
SD309 Ensemble Averager, a Model 13116-2A X-Y Display 
Oscilloscope and an X-Y plotter to provide an instantaneous 
evaluation of the hot wire and tase Sond nee Signals. The 
equipment allowed positive determination that measurements 
were sufficient to define the inertial subrange. A system 
operation diagram is provided in Appendix E along with 
several of the X=-Y plots obtained from the real time 


analysis. 
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ae DATA REDUCTION AND ANALYSIS 


A. PRELIMINARY ANALYSIS 

The analog tapes were played on a strip chart recorder 
to permit initial evaluations of the data. From this, 
data periods were chosen for subsequent reduction and 
analyses. This preliminary analysis disclosed a recording 
problem that had previously gone unnoticed. A periodic 
‘tape transport vibration had occurred due to the 7 1/2 ips 
gear drive and this obliterated the recorded signal. For- 
tunately, there were sufficient time intervals between 
vibrations to permit 20 minute data records to be obtained. 
It was from these data sections that the periods listed in 
Table II were obtained. 

The difficulties encountered in making the temperature 
fluctuation measurements, as described in Section III, 
were evident in the strip chart recordings. The tempera- 
ture data either failed to show sufficient signal strength 
or displayed a DC drift which rendered them impossible to 
digitize. An example of the two cases if§ shown in Figure 
ls 

The paucity of tempernmatune data led to ayfocusing of 
attentilonmomernes=vyelocity Hilluctualtiion m@asur@éments. The 
criteria for selection of subsequent data periods for in- 
clusion in the reduction and analysis phase consisted of: 


Gh) escagnal @eatrengthn (2) averaging time (3) DC drift rate 


oy 








Pago Feet ned equdiwe.Temneratunes Signa ls Showing (A) 
Mislisirelrewemcronaieeatrenyci and (B) DC 
Bie ave t. 


eRe 





and (4) availability of mean profiles. Sufficient signal 
levels were required to minimize amplification during the 
analog-to-digital conversion process. Twenty minute data 
runs proved sufficient to permit ample averaging time for 
Fast Fourier Transform (FFT) calculations. DC drift was 
kept at a minimum to permit ease in digitization and to 
allow greater confidence in the data samples. Finally, 
data periods were selected to correspond to times when mean 
wind information was available. A typical signal meeting 


the above mentioned criteria is shown in Figure 12. 


B. ANALOG-TO-DIGITAL CONVERSION 

Analog-to-digital conversion was performed on a hybrid 
computer system located on the fifth floor of Spanagel 
Hall at NPS. The hybrid system consists of an analog com- 
puter, COMCOR Ci 5000, Hinterfaced electrically with#a 
digital computer, XDS 9300. A detailed discussion of the 
Operatelon of thisPcomputereractlity was given by McKendrick 
(1972) which should be referred to for additional informa- 
E10ne) 8) block didgram, provadedin (igure 13, illustrates 
the analog-to-digital conversion sequence followed in this 


St why. 


C. oy ENS RAC L® NINE~TRACK® CONVERSION 
The analog-to-digital conversion process generated 
seven-track octal base data samples written on digital tape. 


The computer used to analyze the digital data was an IBM 
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Ma gure 15, A-D Conversion Procedures. 
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360/67 computer. This computer utilizes nine-track hexi- 
decimal based digital samples so a conversion between the 
seven-track and nine-track systems was necessary. A dis- 
cussion of theprocedures necessary for making the proper 
conversion to the nine-track hexidecimal based system are 
described in Naval Postgraduate School Computer Center 
Technical Note #0211-08. The two primary programs involved 
in this conversion are the TAPEOUT and CONVERT programs. 
These programs along with a description of digital tape 


operations are discussed in Appendix F. 


D. DIGITAL ANALYSES PROGRAMS 

‘Fast Fourier Transform (FTOR) and spectral analysis 
(SCOR) programs were available in the IBM 360/67 subroutine 
library. The following description of the FFT progran 
package was taken from McKendrick. The programs were 
written by J. F. Garett and J. R. Wilson while students 
at the Institute of Oceanography at the University of 
BrPLecish Columbia. the "OCEAN" “subroutine’s; of the UBCFTOR 
preg tame instructss ches! EMeso0 to read thewinputsGape, one 
block at a time, compute the FFT of the data block and 
then store the resulting Fourier coefficients on the out- 
Put tape. Thescomputer printout flor the ATOR#program is 
amelesting, Cimcacitm@a ug talewohannelimand the Fourier coeffi- 
cients produced. An example of this output is shown in 
Figure 14. 

WitED inese CO@EFIcients, the UBCSCOR program is uséd 


to compute the spectrum for a given channel, or the cross 
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spectrum between a pair of channels, depending on what is 
Selected by the user. The output of the S@OR program con= 
sists of a tabular computer printout, as shown in Figure 

15, and a graphical plot of the spectrum versus frequency 

as shown in Figure 16. Appendix G of this thesis contains 
the preface comments of the FTOR and SCOR programs along 
with additional information for more efficient utilization 
of the spectral analysis package. Several difficulties 

were encountered during utilization of the analyses programs 
and consideration of these difficulties is the primary 


reason for inclusion of Appendix G. 
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Ve RESULTS 


A. ONE POINT SPECTRAL METHOD 

The procedure first selected for estimating € , and 
hence u,, was to use single point values from the spectra, 
Similar to that in Figure 16. A -5/3 slope was approxi- 
mated for the spectrum and the frequency and spectral den- 
Sity, for a spectral value coincident with the line, were 
taken from computer output similar to that illustrated in 
Figure 15. Equation 10 was then used to obtain € from 
values of frequency, spectral density, wimd speed and probe 
height. The computed value of € was then used in Equation 
fe eorObtainea value ier the frictlongwyelocity (u_) eethe 
computer program used for these calculations appears in 
Appendix H. 

OU, Values were obtained by this method) for -all data 
periods and tabulated results for one minute and 10 minute 
averages are listed in Tables III-VII. Fhis method ap- 
peared successful for spectra such as those shown in Figure 
MGpewhieh hak a adéfinite =—5/3 slope for a majority of 
points. However, the single point spectral method is quite 
subjective when used with spectra such as the one shown in 
Figure sl/s eft is eashly seen that ‘selections of unrepre- 
sentative spectral estimates could result in values of € 
Vntemewould Gditfer @s9igqnificantly. Because of the subjec= 


tive manner in which -5/3 slope lines couljdd be drawn and 


47 





TABLE III 


Momentum Flux Calculations for 1103-1138 20 Sept 1973 
uSing the One Point Spectral Method 


Height Time Average € ue 
(m) (min) ane secman (cm sya} 
14.41 1 1.078 S727 
10 as 9 2-129 
T3ne0 1.010 85.32 
10 i359 9.109 
PS, 1 TOIL 9.2.8 741! 
10 iy OMB: 8.865 
4652 1 2.418 7.602 
10 22 O36 7 asl 
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TABLE IV 


Momentum Flux Calculations for 2008-2018 20 Sept 1973 
uSing the One Point Spectral Method 





Height Time Average € u, 
, - be 
(m) (min) (cm sec ah (cm sec 1, 
14.41 AL 9.742 7, O74 
Wt) 8.937 peo 7 2 
8.39 1 6.211 18.291 
10 11.679 15.769 
A 52 ie 20.806 15.643 
10 16.128 4, 287 
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TABLE V 


Momentum Flux Calculations for 2023-2053 20 Sept 1973 
uSing the One Point Spectral Method 


Height Time Average E u, 
(m) (min) (cm? sec™*) (cm sec +) 
14.41 1 9.304 17.505 
10 8.556 #£17.662 
8.39 1 14.042 Nemeg 
10 14.519 16.956 
4.52 1 21.934 15.829 
10 22.634 15.996 


50 





TABLE VI 


Momentum Flux Calculations for 2323-2341 20 Sept 1973 
using the One Point Spectral Method 


Height Time Average atc u, 
(m ) (min) an 5 (cm b yap 
14.41 ali 4.069 14.925 
10 3.343 12.688 ~ 
le rere 1 3.356 12.973 
10 22576 Ll. os2 
8.39 1 7. A65 13.859 
10 6.265 2 .cH3 


on 





TABLE VII 


Momentum er bux Calculations #20r 063/-0652 Al Sept 1973 
uSing the One Point Spectral Method 


Height Time Average € u, 
-(m) (min) vane aaa (cm SA cae 
14.41 1 : 1.923 10.146 

10 e114 8.625 
8.39 1 72330 13.991 

10 5.660 12.386 


BZ 
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Figure 17. Velocity Spectrum Showing Deviation from 
57 oO Lope . 
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a coincident spectral value selected, there was some 


doubt expressed in results from other than ideal spectra. 


B. INTERCEPT METHOD 

To avoid the subjectiveness of estimates of € based 
on single spectral estimates, a second method was utilized. 
This method is based on Equation 10 rewritten in the fol- 


lowing form 


2/3 


Ln kd(k) = -2/3 Ln k + Ln Ge “(14) 


For this expression a plot of Ln kd(k) versus Ln k would 
result in a curve with a slope of one and an intercept 
23 


value equal to Ln ae - From this intercept, the value 


OL € can be obtained from 


Exp b,3/2 

Ss LOR € 
16 J 

where b = intercept value. 


The advantage in using the intercept method versus 
the single point spectral method is that the intercept is 
determined from an average of the points along a line de- 
fined by the -2/3 slope. Figure 18 is an example of most 
curves obtained from the application of Equation 14 to 
spectral@resultsi. Difficulty in approximating the silopé 
was encountered in a few cases such as that illustrated 
ineFigurems.o? but in most instances the f£torcing of a +1 


slope permitted easier detection of spectral harmonics 
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which @Gid not £1 thie =—5/3 critierial. Table XI does not 
indicate a one minute spectral average Since the curves 
for that time average were Similar to Figure 20. It is 
clear that fitting a +1 slope to this scattering of 
points would be difficult. 

The intercept method was used to calculate € and u, 
for all five data periods. The results of these calcula- 


tions are listed in Tables VIII-XII. 


C. COMPARISON OF METHODS 

A comparison of Tables III-VII and Tables VIII-XII 
indicates that somewhat higher u, values were obtained 
with the intercept method. The apparent validity of the 
higher values will be examined later. 

Considering the intercept method, a majority of cases 
vecidedmaargereu. values for longer time averages. The 
importance of resultant higher u, values will be apparent 
when comparisons are made between shipboard results and 
results from other studies. The subjective nature of 
selecting spectral estimates in the one point spectral 
method ~esBmanatestied by larger deviations) in ww. values 
associated with periods when spectra were not "“near-per- 
PeCt tT. 

The one point spectral method was compared with the 
intercept method because of interests in designing pro- 
Gedures Sforequnck jesitimates of wl during measurements. 


dle verluesmor € could be obtained from single spectral 
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TABLEP VILL 


Momentum Flux Calculations for 1103-1138 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average |: a uy, 

(m) (min) (em? Ber cay (cm swe 

14.41 1 1.329 Sy ILieME 

ise So 3). BPE, 

10 e238 o, Ota 

20 1 Pw a 9.248 

eS 12 O 1 1.445 9.139 

4 1.379 8.998 

20 1.503 9.260 

8.39 O35 9.600 

DoS | 9.224 

20 2.622 9.584 

4.52 3.086 8233 

3.148 Sa 2oc 

20 3.689 8.738 
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TABLE IX 


Momentum Flux Calculations for 2008-2018 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average Ee uy 

(m) (min) fom ye (cm per: 

14.41 1 10.288 oO 1 2 

8.285 eel 35 

10 10.935 18.474 

8.39 1 See 26 16.437 

11.828 15.988 

10 13.678 16.622 

4,528 1 i a7 OO 14.731 

14.677 13.846 

10 13 ae eae A 63.5.5 


60 





TABLE X 


Momentum Flux Calculations for 2023-2053 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average € u, 
(m) (min) (ane co) (cm =a 
14.41 1 12.071 19.093 
11.825 iS) aces by 
20 11.540 18.809 
8.39 163723 17.774 
16.496 17.694 
20 6243 17.603 
AZZ 1 25,025 FouD 4 i 
3 24,492 165432 
20 24.149 16.346 
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TABLE XI 


Fad 


Momentum Flux Calculations for 2323-2341 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average € Cis 
(m) (min) (om? ae (cm seo 
14.41 2. 9.2 5 IMS, oye Me! 
10 13.564 19.849 
oe 0 3 36472 To235 
10 14.009 19.488 
8.39 14.418 16.917 
10 5. 025 ee Sal 
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DABLEP AIL 


Momentum Flux Calculations for 0637-0652 21 Sept 1973 
using the Zero Intercept Method 


Height Time Average € Lie 
(m) (min) ae: en (cm fees 

14,41 1 1.209 8.868 
. 1.955 TO. 407 
10 Po570 9.675 
8.39 1 5.35 16 1230 
4 Bo. 451 1 4 SS 
10 y2o43 13.691 
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estimates then the structure function measured may be 
valid. The latter has been shown to be a poor approach 


on the basis of these results. 


D. EVALUATION OF THE VARIATION OF U, WITH HEIGHT 

A long accepted concept in turbulence theory has been 
the existence of a constant-stress layer. This concept 
is often used to evaluate the representativeness of accumu- 
lated data (e.g., Stegen, et al) and will be applied to 
test the validity of assumptions and expressions freed I 
the present study. 

Table I was compiled by Stegen, et al for conditions 
when the measured mean wind speed was 5.12 m 500. This 
mean speed corresponds to the wind conditions for the 
values tabulated in Tables VII and XII. A comparison of 
results in Table I and Tables VII and XII reveals lower 
values of u, in the present study. As previously men- 
EVOnCACgmEme COMSGdaNnt, (Ga) ,.41n Equation LO is omer ica lt 
derived and several values have been suggested in the 
literature. A value of a@ = 0.60 was used in compiling 
the results in Table I while a value of 0.48, suggested 
by Pond, et al (1966), was uSed in the present study. 


U, results listed in Table VIII support a constant- 


* 
flux assumption as well as the validity of Equation 7 for 
Deco lcunietons, £Or nemtyal er near neutral conditions. 


The anomalously low value of u, at 4.52 meters could be 


Enesresult sof acceleration due te the ship"s bow. It is 
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seen that as the spectral averaging time increases the 
value of the flux increases. This would be consistent 
with an interpretation that increased averaging time 
would reduce the observed acceleration effects. 

Tables IX-XI correspond to periods of stable strati- 
fications while Table XII reflects a period of unstable 
stratification. It is seen from Equation 12 that neglect- 
ing diabatic effects would introduce a bias into the data. 
A stability correction could be applied to adjust thie, £lux 


values for the more general case of stratified flow. 


Ee THE VARIATION OF € WITH HEIGHT 

Figur éd 21-25 show the plots of Ln € versus Ln z for 
the periods listed in Tables VIII-XII. On the basis of 
Equation 8, a plot of Ln € vs Ln Z should yield a -1l 
slope for neutral conditions. Results from Stegen, et al, 
Figure 3, were in agreement with this prediction. 

In Figures 22 and 23 the slopes are less than -1 while 
in Figure 24 the slope is greater than -l. An initial 
concern when these results became known was whether possible 
attenuation of the sensors at lower levels was a factor. 
Such a concern is logical since the sensors near the sur- 
face would be likely to lose sensitivity more rapidly as 
a result of increased salt loading at lower levels. 

The data period which occurred last is the one which 
shows the greatest slope which indicates that sensor deteri- 


oration at lower levels was not responsible for the 
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Pigure 21. Measured Values of € vs Height for 1103-1138 
ZC sept 7S FOr Sealing Purposes ja Factom of 
5.0 was Subtracted from the x-Values. 
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Measured Values of € vs Height for 2008-2018 
2ZUmserie i 3. Or Scaling Purpases a* Factor of 
5.0 and 1.4 were Subtracted from the X and Y 
Values, Respectively. 
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Measured Values of € vs Height for 2023-2053 

20 Sept 73. “For Scaling Purposes a Factor of 
9-0 and 2.0 were Subtracted from the X and Y 

Values, Respectively. 
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Figure 24, Measured Values of € vs Heuaght for 0637-0652 
eocp em owt orm —ocalanGg Pwrposes a Factor of 
6.90 was Subtracted from the X Values. 
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Figure 25. Measured Values of € vs Height for 2323-2341 
20) DCD meer Or Scaling Purposes a Factor of 
9.0 and 1.4 were Subtracted from the X and Y 
Values, Respectively. 
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variation in slopes. Further support of this conclusion 
is provided by using Equation 11 as an imdependent method 
of determining the dissipation. 

For the 2023-2053 20 Sept 73 data period, the spectra 
at all levels were similar to that shown in Figure 26, 
where scale separation between the pure inertial transfer 
and the viscous dissipation meq teres was discernable. The 
upper limits of the inertial subrange were determined from 
the spectra and Equation 11 was used to obtain an estimate 
of the dissipation. In Table XIII the values obtained by 
this method are compared to those obtained from spectral 
estimates. The inner scale computation iS independent of 
sensor calibration Since it is determined by frequency dis- 
tributions and not on the spectral value. The apparent 
agreement between the two estimates is am in situ calibra- 
tion which indicates the loss of sensitiwity by individual 
sensors was not the cause for the deviation from the pre- 
dicted slope. 

The deviation from the predicted slope was also ex- 
amined for possible influence of non-neutral stratifica- 
tions or wind-wave coupling. Deacon (1951), from examina- 
tions of numerous vertical wind profiles, proposed the 


following expression for the wind shear 


Ju uy, Z, = 
+t = | B (15) 


where B is a function of stability. 
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Figure 26. Twenty Minute Velocity SPCOGerumerors 2023-2053 
20 Sept 73. 
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TABLED XE Eo: 


A Comparison of Dissipation Obtained by Spectra (€) and 
Dissipation Obtained from the Inner S carlic (e) 


Height E E U, 
(m) - 2 ~3 2 = 


14.41 Slope oe DG 17.662 
a3 9 14.519 14.304 16, 736 
4.52 22.634 22.492 i ..996 
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The variation of $8 as a function or Ri, proposed by Dea- 
con, is shown in Figure 27. In Deacon's results, the 
variation of 8B may be described in terms of hydrostatic 
stability. For neutral conditions, $8 = 1 so Equation 15 
then reduces to Equation 6 and the -1l slope on the Ln € 
versus Ln z plot is predicted. 

For unstable conditions §$ is greater than 1 and for 
stable conditions § is less than 1. Deacon's results 
predict that the expected slope of a Ln € versus Ln 2 
plot would deviate from the neutral case with slopes less 
than -l for periods of stable stratifications and with 
slopes greater than -1 for unstable conditions. Examina- 
tions of Figures 22-24 indicate why such considerations 
are important in interpreting the shipboard results. 
Figures 22-23 correspond to periods of stable stratifica- 
tions with the 2023-2053 20 September data period being 
the most stable. In both figures, slopes less than -l 
occur, with the more stable of the two periods showing 
the greatest deviation — the -l slope. The results 
are, therefore, consistent with arguments proposed by 
Dieacem, Figure 24 corresponds to an unstable period and 
shows a slope greater than -l as predicted by Equation 15. 

The results were further examined with respect to pre- 
dictions from the Monin-Obuhkov similarity theory. The 
theory states that in the boundary layer there exists a 
characteristic length determined by the stability and tur- 


bulence which can be used as a scaling parameter. OL 
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Figure 27. Variation of 8 with Ri (after Haltiner and 
Martin). 
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neutral conditions the appropriate characteristic length 
is the height above the marine boundary. 

Equations 12 and 13 show the effect of non-neutral ; 
considerations in developing an expression eee eins the 
vertical variation of €. Equation 12 indicates a relation- 
ship between the scaling parameters Ri amd L. Such a 
relationship was described by Busch (1973) ina review of 
information available in the literature. As presented by 
Busciwetine Lunction a) 5 in Hquation 12 has’ thle form 
Ge = =) where several flux-profile relationships have 
been suggested for >a Table XIV includes the formula- 
tions which Busch used to describe ad along with constants 
estimated by Paulson (1970). Relationships between the 
Monin-Obuhkov length and the Richardson mumber were first 
used to obtain stability length (L) corrections for the 
shipboard results, since Richardson numbers were available 
from profile meaSurements described by Cavanaugh (1974). 
Values of L were then used to define (>. = = from rela- 
tionships in Table XIV. 

Figures 28 anc 29 show the effects of stability correc- 
PLoS Seater sper iods with StabMes Stratgedcations listed 
in Tables IX and X, respectively. Comparisons of the cor- 
rected andwuncormected Tesults indicate) ¢hat adjustments 
were toward the predicted -l1 slope. The fact that the 
Stability correction produces a -l1 slope supports the basic 
premise of the Monin-Obuhkov theory and the validity of 


the welationships in Table XIV. 
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TABLE XIV 


Flux-Profile Relationships Describing the Stability 
Function, on? 


9 a Stability Condition 
Z Z 
+7— — > 
ae) = - > 0 
Zz =1/4 7, a 
-~16= — < 
(1-165) = 0 
(isi ari) 17" ao SOs Il 
CE Sy -O.1 < Ri < 0,036 
0.88(1-6Ri) 0.036 < Ri < 0.1 
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Figure 28. Effect of Stability Correction for 2008-2018 
2Ue Sept 73. 
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Figure 29. Effect of Stability Correction for 2023-2053 
ZO Sept. 73¢ 
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The Monin-Obuhkov theory and the flux-profile relation- 
ships were also examined for the period with unstable stra- 
tification. Figure 30 shows the results of this examination. 
Poo) See iee-in Oa ro thewesis sa tendeney for#adjustment 
toward the predicted -1 slope. However, further decreases 
in the Richardson number caused an unexpected reversal of 
the adjustment. 

It is seen from Figure 30 that a value of -1.0 for the 
Richardson number results in a vertical variation of € which 
deviates significantly from the predicted slope. The reason 
for thits’ deviation can BbeWexplained by examining the func- 
don, ()_ (— - =). Figure 31, from Garratt (1972), shows a 
plot of this function uSing the expressions in Table XIV. 

To provide the correction necessary to obtain a -1 slope, 
(9, - =) must decrease for decreasing values of a It is 
apparent from Figure 31 that this criterion is not met. On 
the basis of these results there 1S some evidence that the 


1/4 


expression p  (1-16=) cannot provide the adjustment 
necesSary tc correct the plot for conditions of large 
instability. 

The fact that the: stability corrections did not com- 
pletely restore the -1 slope, even in the stable case, 
leads to examinations of the possible effects of wind-wave 
Coup! ino Davidson and Frank (1971) observed wave-related 


velocity fluctuations and wave-related momentum transfer 


in spectral results obtained over natural waves. Davidson 
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Figure 30. Effect of Stability Correction for 0637-0652 
2A Sept 7 3. 
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Figure 31. Non-dimensional Dissipation Rate as a 
Bunet1ton Of 2/L G@atter Garratt) . 
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(1974) described the influence of wind-wave coupling as 
having the same effect as would be due to the influence 
of stable stratification. The wave-related effect de- 
scribed by Davidson was a reduction in momentum transfer 
from that associated with neutral conditions. The re- 
sults of the present study are consistent with Davidson's 
findings in that there is a decreased momentum transfer, 


u as the surface is approached. If this were due to the 


x? 
influence of wind-wave coupling, it would yield the ob- 


served deviation from the neutral slope, shown in Figures 


22and 23. 


F VARIATION OF TEMPERATURE FLUCTUATIONS 

Figure 32 represents the only temperature spectrum 
Obtained from the five data periods which displayed the 
predicted -5/3 slope. Figure 33 is an example of the 
temperature spectra obtained from the majority of data 
runs. The paucity of temperature data precluded serious 


evaluations of the temperature fluctuation expressions. 
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Figure 32. Ten Minute Temperature Spectrum for 2323-2341 
20 Sept 73. 
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Figure 33. Ten Minute Temperature Spectrum for 2008-2018 
ZoeSepe 73. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The results obtained from the one poimit spectral 
method are generally lower than those obtained by the 
intercept method. However, the simplicity of the one 
point spectral method would permit a more rapid deter- 
mination of the dissipation and if a general estimate 
were all that was desired, this method would be accept- 
able. The results of the intercept method show excel- 
lent agreement among levels and with results from other 
studies. The longer time averages (~20 min.) appear 
to be necessary to avoid acceleration effects at the 
bower kevel. 

The close agreement between dissipation values de- 
rived from the spectral estimates and the inner scale 
calculations confirm the spectral estimates and also the 
enoilce sof the constant, a = 0.48. Results on the varia- 
tion of € with height also verified the formulations 
used in calculating the dissipation and flux. In addi- 
tion, the sensitivity of the results to the varying con- 
ditions of hydrostatic stability and, perhaps, wind-wave 
coupling support the predictions proposed by Deacon and 


Davidson, 
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B. RECOMMENDATIONS FOR FUTURE EXPERIMENTS 

It is strongly recommended that a real time spectrum 
analyzer be used in future observational experiments. In 
addition to the capability for instantaneous evaluation 
of spectral characteristics, the resulting spectral plots 
could be used for general estimates of the dissipation. 
If observational experiments are to be made on a more fre- 
quent basis, there is insufficient time to properly 
evaluate the data before the next experiment. The use of 
real time analysis techniques would preclude the need for 
digitization and tape Manipulations and save valuable and 


costly computer time. 
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APPENDIX B 


A. CALIBRATION AND SCALING PROCEDURES 
1. Velocity Calibration 

Calibration of the velocity sensors was accom- 
plished using a manometer and a TSI Model 1125 Calibrator. 
The calibration procedure was to record woltage output ver- 
SuS manometer height over a certain range of velocities. 
The computer program included in this Appendix was Berti sed 
to generate the calibration curve for further uSe in estab- 
lishing the calibration factor. This procedure was followed 
for each of the velocity systems and resulted in the cali- 
bration curves included in this Appendix. 

The calibration curve obtained from this program 
resulted ina plot of v? versus Yu where V is voltage and 


u is wind Speed. 





LE Intercept = b 


90 





This curve describes a line 


2 


V= arcana <a 


oD (1) 


For use in data analysis it was necessary to obtain a 


scale such that 


c xe’ =u' 
where C = calibration factor; e'.= voltage fluctuation 
and u' = velocity fluctuation. Differentiation of Equa- 
tion (1) yields 7 
ey Aen) Oa a 
or 
1/2 
SE WU i eee Oe 
a 
or 
1/2 
pees ar eet 
a 


The calibration factor is then 


4vur/? 


C a cel 


where "C" scales between volts and cm/sec. 

The procedures for actual calculation of this 
factor consisted of (1) obtaining the slope of the cali- 
bration curve (2) selecting a "u" valve equal to that 
observed during the measurement period and (3) entering 
the curve with this "u" value to obtain the corresponding 


voltage. 


bed 





2. Scaling Procedures for Estimation of Epsilon 


In terms of the output of the spectral analyses 
programs the units associated with f$(f) or kod(k) were 
en st) ce An examination of Equation (1), Section II, 
shows that the units should be FES - Recalling that 
the calibration factor previously discussed scales between 


volts and cm/sec, it is necessary to multiply the € value 


by 
2 2 3/2 3 
(cm /sec ) cm/sec 3 
ead = a at? (C) 
volts 


eee 1 Cees 


to obtain the value of € in ery ae can 

To obtain the correct value of € it is also neces-~ 
sary to consider the gains used during the record and 
analog-to-digital conversion phases. In addition a scale 
factor is needed that scales between the analog voltage 
and the digital output. A 10db gain was used when the 
raw signals are recorded and an additional 50 db gain was 
employed during the A-D process. The scaling between the 
voltage input and the corresponding digital output is ob- 
tained by considering that the maximum voltage input to 
the COMCOR Ci 5000 is 100 volts and this corresponds to 
a digital maximum of 273 on the XDX 9300. 


The total gain and scaling for the observational 


data then becomes 


digital Value obtained 
UO) GL RCE Rar ae from IBM 360 
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where V = signal from velocity unit 
Gl = amplification prior to record (10 db) 
G2 = amplification during A-D comversion (50 db) 
SC = analog-to-digital scaling factor (273/100) 
To obtain the value which corresponds to that which can 


be used with the calibration factor previously obtained 
V = (IBM 360 value x 1/SC x 1/Gk x 1/G2) 
The final scaling then becomes 
3 
jiscaie x Gainl x Gain2] x c | 
or in this instance 


2 : 
t100/2 een x Of O02] x c | 
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Figure El. Real Time Spectral Analysis of Velocity Signal 
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Figure E2. Real Time Spectral Analysis of Velocity Signal 
EOm~e2O023—-2057820 Septm73. 


103 








APPENDIX F 


A. TAPEOUT PROGRAM 

TAPEOUT is a general purpose program used to determine 
the label and file contents of magnetic tapes. The program 
can take an otherwise unreadable dale and provide the user 
with Rares cient information to permit access to the tape. 
TAPEOUT will Prints a £1gle number, density, parity, length 
in bytes of the record and a dump of the record. In this 
study the two key outputs were the record length and record 
dump. These outputs permitted the user to obtain the block- 
size and ascertain the location of erroneous records. The 
following is an example of the TAPEOUT deck set-up used in 
this study. 

//JOHNSTON JOB (2413,0392,XM24), ‘JOHNSTON SMC 2561' 

//EXEC TAPEOUT, PARM='2.0,,2' . 

//TAPEIN DD UNIT=2400-1,LABEL=(,BLP), VOL=SER=KARLYN 

/* 

TAPEOUT is part of the IBM 360/67 subroutine library 
and a discussion of the parameters (PARM=) is available in 
the computer center consultant's office and also in Com- 


puter Center Technical Note No. 0211-08. 


B. CONVERT PROGRAM 
The seven-track octal tape is converted into nine-track 
hexidecimal data samples usince the CONVERT program. The 


DEQCKEGCHaLaCeeristics 6f@the data samples are still 
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Maintained on the nine-track tape but two additional para- 
Meters are affixed to each block. These two parameters, 
NCHAN and KMAX, are affixed to the beginming of each block. 
NCHAN is the number of channels of analog data digitized 
on the seven-track tape. KMAX is a numerical value set 
equal to the maximum number of samples per block. The in- 
Ssertion of these numerical parameters at the beginning of 
each block and the change in the block format by the FORM 
Subroutine will put the seven-track tape in the proper for- 
mat for use in the UBCFTOR program. The CONVERT program 
used in this study along with sample JCL cards is given 


below. 


//TOHNSTON JOB (2413,0392,XM24),'JOHNSTON SMC 2561' 
J//CONVERT EXEC FORTCLG,REGION.GO=100K 
//FORT.SYSIN DD 
DIMENSION IDAT(1024), DAT(1024) 
FACTOR =* 1), 
REWIND 2 
REWIND 4 
NCHAN 
NRECL 
J = 0 
10 READ (2,20,END=60,ERR=40) IDAT 
20 FORMAT (20(50A4) ,24A4) 
J=J+1 
CALL FORM(IDAT,NRECL) 
DO 30 I=1,NRECL 
WRITE (4) NRECL,NCHAN, DAT 
GO TO 10 
40 WRITE (6,50) J 
5O FORMAT ('0' ,5X ,)READJERROR, RECORD NO.=',14) 
GO TO 10 
60 WRITE (6,70) J 
(Om FORMaAwwe (9G) ok, END SOF TAPER REGORD NO.=' 714) 
END FILE 4 
RETURN 
END 


& 
1024 


/* 
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//GO.FTO2FOO1 DD UNIT=2400-1, VOL=SER=KARLYN,LABEL=(,NL), 

tad, DISP=OLD, DCB= (DEN=2, RECFM=VS, BLKSIZE=8192) 

//GO.FTO4FO01 DD UNIT=2400, VOL=SER=NPS170, DSNAME=DEREK, 

LABEL=(1,SL) ,DISP=(NEW,KEEP) , DCB=(DEN=2,RECFM=VS, 

Me A BLKSIZE=8200) 

It is noted that the user inputs to the program consist 
of the dimension of the arrays IDAT and DAT and the value 
of NRECL. IDAT and DAT should be dimensioned at least NSAMP 
x NCHAN thus signifying the need for recording the digitiz- 
ing parameters. The argument NRECL must be assigned a 
value equal to NSAMP x NCHAN. The value’ of the scaling fac- 
tor has been set equal to one in this study since it was 
found to be easier to perform all calibrations and scalings 
in one step following the use of the spectral analysis pro- 


grams. Computer Center Technical Note No. 0211-08 provides 


additional information concerning the CONVERT program. 


oe DIGITAL TAPE OPERATIONS 

Unlike the Hybrid computer system, where a].l tape mani- 
pulations were performed by the user with technical assist- 
ance from the computer FaASits staff, the IBM 360 computer 
system manipulations were performed by the staff members. 
In order to have the tapes handled as desired, explicit in- 
structions had to be given. These instructions were given 
in two ways. The first instruction consisted of a Job Re- 
quest Form addressed to the computer center staff identify- 
ing the tapes to be utilized and whether they were to be 


read off of (RINGOUT) or written on (RINGIN). The second 
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instructional procedure was conducted usimg job control 
language (JCL) addressed to the computer and read as part 
of the submitted program. 

1. Job Control Language 

Job control language cards provide information to 
the computer about tape identification, tape disposition, 
tape data files, tape mount anaes ae and data format. The 
cards are appended to the program and read through the “hot 
card reader" as part of the program, whereas the Job Request 
Form is submitted "over the counter" to the computer center 
personnel. An example of a JCL card for ttape processing 
would be: 

//GO.FTO8FO0O1 DD UNIT=2400-1,VOL=SER=KARLYN, LABEL=(1,SL), 

VAIS DSNAME=DEREK,DISP=(NEW,KEEP), DCRB=(DEN=2,RECFM=VS, 

// BLKSIZE=8192) 

While the majority of the JCL Pee remaims unchanged, there 
are portions which are subject to change trom tape to tape. 
These portions are discussed below. 

FTO8 - This group indicates the tape mount on which the 
tape is to be placed. The two digit number (08 in this case) 
must match the number provided in the pragram's "READ" or 
"WRITE" statement. 

FOO1l - This group indicates the sequential number of 
passes Enwough thie tape which has occurge@ up to this point. 
In this case it is the first pass. 

UNIT=2400-1 - This group specifies the tape as a seven- 


track tape. The nine-track designation ais 2400; 
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VOL=SER=KARLYN - This group identifies the tape to 
be processed. In this case the name "KARLYN" signifies 
the tape is an "external" tape with the mame "KARLYN" 
Physically appearing on the tape. An “imternal" tape 
would have a serial number such as NPS211. These tapes 
are kept on file in the computer center and cannot be 
Physically removed from the facility. 

LABEL=(1,SL) - This group identifies the tape file to 
be processed and specifies that the tape has a standard 
label. "External" tapes may be processed as no label 
tapes (NL) but “internal" tapes will be given a standard 
label before the user gains access to them. This is accom- 
plished by the computer center staff processing a "canned" 
program which places a standard label on the serialized 
tape. 

DSNAME=DEREK - This group specifies the data set name 
(DSNAME) applicable to the specified LABEL group. The 
data set name may be a unigue name limited to six alpha- 
numeric characters. Once a file has beem given a data set 
name access to the data file can only be obtained if this 
same DSNAME is used. This process can be circumvented by 
the use of by-pass label processing (BLP) but does provide 
protection against accidental erasures. 

DISP=(NEW,KEEP) - This group specifies the tape dis- 
position once it has been processed. 

DCB= (DEN=2, RECFM=VS, BLKSIZE=8192) - This group is re- 


ferred to as the data control block. The density (DEN=2), 
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tape record format (RECFM=VS) and block size (BLKSIZE= 
8192) are stipulated uSing this group. The density speci- 
fies the number of bytes per inch, in this case the 

number 2 denotes 556 bytes per inch. This corresponds to 
a manual setting of the DENSITY dial on the face of the 
seven-track tape unit. RECFM=VS specifies that the record 
format is variable spaced, the result of recording and 
digitizing on the seven track unit. The block size speci- 
fies the number of bytes per block of data. If 2048 words 


per block were used, four bytes per word would result in 


8192 bytes per block. 
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APPENDIX G 


A. GENERAL 

All the UBC time series analysis programs (FTOR, SCOR, 
FC PLOT) were originally listed on faculty tape NPS216, 
file 1-3 respectively. The Bilig nadie have since been placed 
on computer center disk, MARY, and are stored in machine 
language which disallows a ii Eine of the program for user 
Study. However, source decks for UBCFTOR and UBCSCOR have 
been retained by Dr. Kenneth Davidson of the Meteorology 
Department anc by Dr. Noel Boston of the Oceanography De- 
partment. These source decks can be used to obtain a 
listing of the programs for individual study. The most 
recent computer center documentation on the use of the time 
series analysis package is given in the following memoran- 
aumm@wiittitien by Miss) Sharon Raney of the computer center 
staff. This memorandum provides the information and JCL 
necessary for utilization of the disk-stored analysis 


programs. 


B. Er LOR 

The preamble to the FTOR program is given following 
this discussion. It provides the user with the informa- 
tion needed to prepare the input data cards necessary for 
program execution. One of the most misleading parameters 


called for in the input values is the variable NBLOCK. 
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This value is provided by the user on the second data 
card with the value supposedly entered in columns 4-5. 
In actuality the value is read as a five digit integer 
thereby making the maximum value of NBLOCK 99999 vice 
99 as implied by the program preface. | 

An important consideration involves the value of the 
variable MPRINT. This variable controls the computer 
printout of the FTOR program. If the user desires to sup- 
press the printout the insertion of a specific value can 
save considerable paper and cut computer time Mert cancion 
In the unsuppressed mode, the output will follow the format 
shown in Figure 13, Section IV. The suppressed output will 
consist of one line outputs for each data block which will 
Mimo at iemiccr that scoctficientsstor block ~ have been 
Vea Me CellC SeOt atic souLDUutL Capen a Bysusings the 
suppressed output, computer turn-around time was signifi- 


cantly reduced. 


C. = G0). 

The preamble to the SCOR program also provides the user 
with the information necessary to compile the data input 
cards. The SCOR program requires the input of several 
variables whose values are determined during the A-D con- 
version process. One of the values on the first data card 
is the number of data blocks on which the spectral analysis 
will be performed. This value is set by the user and can- 


not exceed the number of data blocks processed by the FTOR 
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program. As was the caSe in the FTOR program the maximum 
number of blocks which can be handled is 99999 and not 99 
as the preface would lead the user to believe. 

As was discussed earlier in this thesis the output of 
the SCOR program consists of a computer printout and a 
graphic plot. The graphic plot consists of an output in 
which spectral density is plotted against frequency for 
all channels being processed. Since the original formula- 
tion of the SCOR program, the computer facility at the 
Naval a ee. School has installed a new plotting 
routine called NEWPLOT. This routine limits the usSer to 
30 inches of plotting paper. The author found it necessary 
to submit a job request asking that plotting routine CAL- 
PLOT be used to prevent early termination of spectral 
plotts@eonsisting of four or more channels. The’ computer 
center staff has been alerted to this problem and any 
future user should consult with a staff member before 


using the SCOR program. 
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8 June 1973 


MEMORANDUM 


From: Sharon D. Raney 
To: Users of UBC Time Series Analysis Package 


Subj: New Operating Instructions for UBC Time Series Analysis Package 
1. The new FTOR deck setup is given in Figure l. 


The dimension size of array DATA must correspond to the block size 
of the data on the input tape. This tape is generally the data as 
converted from the seven track XDS generated tape to the nine track 
IBM/360 tape. 


Formerly, this array was dimensioned in subroutine OCEAN2 and had 
to be changed whenever the sample size changed. This meant the 
entire source deck for FTOR had to be used for each run. 


FTOR is now compiled in FORTRAN IV (H) which generates a more 
efficient code and consequently speeds up execution time. 


2. When the same dimension size is to be used repeatedly, Figure 2 
shows how to obtain an object deck and how to execute FTOR with the 
object deck. 


3. Figure 3 shows the deck setup for SCOR. 


4. Notice the change in the format of the Job Control Language (JCL) 
to define the tape files. These changes enable more efficient 
processing by the operator. Note the use of the IN parameter in 
the LABEL parameter. This must be specified when the file is to 
be read only. OUT is used when the file is to be written only. 
All data sets are PASSed not KEEP. The UNIT-AFF is used on sub- 
sequent files of a tape as is the VOL=REF. 


5. Any questions concerning these new operating instructions should 
be directed to Sharon Raney at X2631 or Ingersoll Hall Room 107. 
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